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Abstract 

We analyze in detail lepton-flavor violation (LFV) in the charged- 
lepton sector such as /i — > 67, r — > /i7, ^ — > eee and the /i — > e 
conversion in nuclei, within the framework of supersymmetric models 
with lopsided Froggatt-Nielsen charges, in which the large mixing in 
the neutrino sector as well as small mixings in the quark sector can 
be naturally accommodated. We show that the present experimental 
limits on the LFV processes already exclude some of the models. The 
future proposed search for LFV, especially in muon processes, can 
provide a significant probe to this framework. 

We also stress the importance of the measurement of U^l^^ in 
neutrino experiments, and the fact that the KamLAND experiment 
could play a significant role to test a certain class of models. 



1 Introduction 



One of big mysteries in the Standard Model (SM) of elementary particles 
is a problem of fermion masses. Since Yukawa couplings, which determine 
the magnitude of fermion masses, are totally free parameters in the SM, 
theoretically we do not know how we can predict a wide variety of masses of 
fermions. 

On the experimental side, evidence of non-zero neutrino masses from the 
atmospheric neutrino experiment has recently been announced by the Super- 
Kamiokande collaboration [|l| . This result is very interesting because not only 
it suggests non-zero neutrino masses, but also it indicates a large mixing in 
the neutrino sector. The tiny but non-zero neutrino masses clearly imply new 
physics beyond the SM, and the large mixing in the neutrino sector suggests 
that a flavor structure in the lepton sector seems to be very different from 
that in the quark sector. Therefore, finding the unified picture between small 
mixings in the quark sector and large mixing in the lepton sector will be an 
important key to understanding the problem of fermion masses, and a lot of 
attempts have been made 0, |^. 

One of the interesting and simple mechanisms to realize the mass hier- 
archy of fermions is the Froggatt-Nielsen (FN) mechanism [Q, which uses a 
broken U(l) family symmetry. It has been proposed that lopsided FN U(l) 
charges for lepton doublets would be an interesting candidate to naturally ac- 
count for the large mixing for neutrinos and the small mixings for quarks . 
The low-energy consequence in neutrino physics has been studied in Ref . . 
It has been also interestingly pointed out that this framework can explain a 
baryon asymmetry in the present universe [Q. 

The supersymmetric (SUSY) extension of this framework, assuming SUSY 
is broken at a high-energy scale, is rather interesting, because we can expect 
the other low-energy consequence. The lopsided structure of lepton doublets 
induces a large mixing in the Yukawa matrices of leptons. Such a large 
mixing has a potential to generate the large lepton-flavor violation (LFV) in 
slepton masses through the renormalization group (RG) effects ||^. Then 
LFV processes in the charged-lepton sector such as /i — *• 67, r /i7, 
/i eee, and the /z — > e conversion in nuclei, are induced through diagrams 
mediated by the sleptons. In the presence of the large neutrino Yukawa 
coupling, the event rates can be within the reach of future experiments or 
the models can be strongly constrained [^, H, H, p!0| . 

In this paper, we analyze the LFV in SUSY models with lopsided FN 
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U(l) charges in detail. We show that the search for LFV, especially in muon 
processes, provides a great impact on this framework, and even at present 
many of the models are almost excluded. The future proposed experimental 
improvement of the search for LFV in muon processes will be significant for 
the SUSY models with lopsided family structure. Therefore, we emphasize 
that the LFV search would be an important window to seek for the answer 
to the fermion mass problem. In section 2, we briefly introduce the models 
with lopsided FN U(l) charges. In section 3, we discuss neutrino masses 
and mixings in two interesting classes of models, and especially stress the 
importance of the measurement of the neutrino mixing U^^^ and the fact 
that the KamLAND experiment could play a significant role to test a certain 
class of models. In section 4, we discuss LFV in detail and show that the 
/X — *■ 67 process is more sensitive to the models than r — >■ /i7. Present and 
future experimental searches have a great potential to probe the models. 

2 Models with lopsided Froggatt-Nielsen U(l) 

It has been pointed out that lopsided FN U(l) charges for left-handed lepton 
doublets Li {i = 1-3) are interesting possibilities to explain the large u^- 
Vr mixing observed by the atmospheric neutrino experiments as well as the 
mass hierarchy of charged leptons and quarks. Here we briefly introduce two 
interesting classes of models. 

In order to account for the tiny neutrino masses, we consider the seesaw 
mechanism introducing heavy right-handed neutrinos Ni {i = 1-3) |TlI|. One 
possible interesting U(l) charge assignment is that the lepton doublets of 
the second and third families {L2 and L3) have the same U(l) charges r and 
the first family Li has a different U(l) charge t + 1, while the right-handed 
charged leptons Ei {i = 1-3) have U(l) charges 2, 1, 0, respectively. We refer 
to this class of models as "Model I" . Another interesting charge assignment 
is that all lepton doublets Li have the same charge r and the right-handed 
charged leptons Ei have charges 3, 1, respectively. Below, we refer to this 
class of models as "Model 11". We list the FN U(l) charges of Models I and 
II in Table |l|. 

Mass terms for the lepton sector are given by 

W = EJ'^LjH, + NJl^LjH^ + ^N,M,,N,. (1) 
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After diagonalizing the charged-lepton Yukawa matrix (/g) and right-handed 
neutrino mass matrix (M) and taking into account the FN charges, we obtain 
the following mass matrices: 



me. = 



feV cos/3 
V2 



diag.(eie-^ 626, 63) e^ma, 



M = diag.(nie^^ n2e'^\ e^^^) Mr, 



fvV sin j3 



(2) 
(3) 

(4) 



where v = ^ [HuY + {H^Y, tan/3 = (Hy) /{Ho), and ms and Mr represent 
a weak scale and a right-handed neutrino scale, respectively. The coefficients 
Cj, Aj, Bi, and are undetermined but expected to be of order 1, and 
5 = and 1 for Models I and II, respectively. Through the seesaw mechanism, 
assuming Mji ^ TO3, tiny neutrino masses can be obtained: 



Cfe^' B,ae' AQe' 
Bide^ Bf AB 
AiCiC^ AjBi 



m 



3 ,2t 



M, 



R 



2 



A 



where coefficients Ai, Bi, and Q are given by 



(5) 
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Table 1: Proggatt-Nielsen charges for matter and Higgs fields. In SU(5) 
language, lOj = {Qi, Ui, Ei), 5^ = {Di, Li), and H is for all Higgses. 
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As can be seen in Eq. (^) , a large mixing between z/^ and Vr can be expected 
since the matrix elements {mi,)ij = 2,3) are of the same order. Note 
that both Models I and II have the same hierarchical structure in charged- 
lepton masses (Eq. (|])); on the other hand, the neutrino mass matrix Eq. 
depends on the models (5). It has been studied in Ref. P that in order to 
obtain the correct masses for charged leptons, the best value for e is 0.07. 
Therefore in our analysis, we will fix e to 0.07. 

The neutrino mass matrix in Eq. (^) is diagonalized by a mixing matrix 

jjMNS. 

[/^^^^^m,?7*'^^^ = diag.(m,i, m,2, m.a), (7) 
^Fa = U^'uM^, (8) 

where and um are the flavor and mass eigenstates of neutrinos. The Dirac 
neutrino mass matrix Eq. (^) is diagonalized by the bi-unitary transforma- 
tion: 

F^^^'^^^m.z^f/''''"^ = diag.(m,Bi, m,D2, m.^a)- (9) 

The mixing matrix f/-^*'"'*^ is relevant to LEV in the charged lepton sector as 
we will see later. Although in general the matrix JJ^"'""^ is different from the 
neutrino mixing matrix f/A'-f^^s 

an important point is that the large mixing 
originates from the lopsided structure in Dirac neutrino masses, and hence 
both mixing matrices JJ^^^^ and JJ^"'"''^ possess a large mixing. This fact is 
very important for the LEV phenomena. 

We should stress that because of the lopsided EN charges listed in Table |l|, 
the Kobayashi-Maskawa matrix in the quark sector has only small mixings. 
Therefore, this framework can naturally accommodate small mixings in the 
quark sector as well as large mixing in the neutrino sector. 

3 Neutrino masses and mixings 

As can be seen in Eq. (^), the mass matrix for neutrinos does not depend 
on the EN charges of right-handed neutrinos. A dependence on the EN 
charge r is found only in the overall factor of the neutrino mass matrix 
Eq. (^). Therefore the predictions for the ratio of neutrino masses and the 
neutrino mixings are almost the same in the different choices of the charge 
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T However, some of the matrix elements depend on 5. In this section, 
therefore, we consider Models I and II separately. 



3.1 Model I {5=1) 

First let us consider Model I. When we diagonalize the matrix in Eq. 
with 5 = 1, we naively obtain the neutrino mixing matrix JJ^^^^ as follows: 

/ 1 0(e) 0(e) 

0(e) 0(1) 0(1) I, (10) 
\0{e) 0(1) 0(1) 



U 



MNS 

r\-/ 



in the leading order of e. Since the matrix elements {m^)ij = 2,3) are 
of order 1, we can naturally get a large mixing for atmospheric neutrinos. 
For solar neutrinos, it is easy to get a small mixing since e2 can be 

naturally small. However, it is not difficult to get even the large mixing 
solutions for solar neutrinos as we will see in our numerical analysis later. 
An interesting point is that a component U^^^^ is expected to be of order 
e, so that this is not only consistent with CHOOZ and atmospheric neutrino 
experiments, but also can be within the reach of future neutrino experiments. 

In our numerical analysis, we randomly generate data sets of coefficients 
y4j, and Oj, in which we vary the absolute values of the coefficients Ai, 
Bi, and Ci from 0.5 to 2 and their phases from to 27r since the coefficients 
Ai, Bi, and Ci are complex constants of order 1. For simplicity, we fix rtj 
(i = 1, 2) to 1 in right-handed neutrino masses in Eq. (^. Then we calculate 
neutrino masses and mixings, and look for neutrino solutions that will solve 
atmospheric and solar neutrino problems. In our analysis, we impose the 
following conditions: 

sin^ 2^,tm > 0.8, (11) 
for the atmospheric neutrino solution, 

0.2 < tan^^soi < 1, 

3.3 X 10-=^ < < 10-\ (12) 



^ Running effects in the RG equations depend on the FN charges of right-handed 
neutrinos and left-handed leptons. Thus the neutrino mass matrix at the low-energy scale 
can be affected by the runnings from GUT scale to right-handed neutrino mass scale. 
However, we have checked that these effects are not very significant, so the results for 
neutrino masses and mixings are almost the same in the different cases. 
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for the large mixing angle (LMA) MSW solution, 



2 X 10"^ < tan^^soi < 3 x 10"^ 



1.3 X 10"'' < — ^ < 3.3 X 10-^ (13) 



for the small mixing angle (SMA) MSW solution, 

0.3 < tan^ ^soi < 3, 



1.3 X 10"*^ < — ^ < 1 X 10~^ (14) 



atm 



for the LOW or vacuum solution. Here Sml^^ = ml^ — '^"^soi = 

m^2 ~ ^liy ciiid we define "effective mixings" sin^ 26'atm and tan^6'soi as 

sin^2^^.,^ = 4|f/,r"P(l-|f/.r"r), (15) 
tan^^^.,, ^ \Ur'/Ur'\^ (16) 
for three-flavor neutrino oscillation. We also impose the CHOOZ constraint: 

\U^^^^\ < 0.15. (17) 

To fix a scale for neutrino masses, we take the scale for atmospheric neutrinos 

as 

5ml^ = 3 X 10-3 eV^ (18) 

for simplicity. 

With these constraints, 0.1% of data sets out of all data sets we generated 
passed the constraints for the CHOOZ, atmospheric and SMA solution, 
0.3% for the CHOOZ, atmospheric and LMA solution, and 0.07% for the 
CHOOZ, atmospheric and LOW solution. Therefore this class of models can 
accommodate various solar neutrino solutions (LMA, SMA and LOW), but 
not the vacuum solution. In Fig. |I] we show the distribution of the neutrino 
solutions. All points satisfy the atmospheric neutrino constraint in Eq. (pT]), 
as shown in Fig. |I[ The circle points satisfy the condition for the LMA 
solution in Eq. (|12|), the diamond-shaped points for the SMA solution in 
Eq. (0), and the square points for the LOW solution in Eq. (|T3|). 

One interesting point is that most of the predicted values for \U^^^^\ 
are larger than 10~^ (Fig. ||), so that the values \U^^^^\ can be reached by 
future neutrino experiments |T^. Therefore the future precise measurement 
of \U^^^^\ will be very important to test this class of models. 
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Figure 1: Distribution of the points of neutrino solutions in Model L Black 
circles satisfy the constraints for the LMA solution, the diamond-shaped 
points for the SMA solution, and the square points for the LOW solution. 
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Figure 2: Predicted values of \U^l^^\ in Model I 
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3.2 Model II {5 = 0) 

In this case, all matrix elements in Eq. (^) with 6 = are of order 1, and 
hence we can expect all elements of the neutrino mixing matrix to be of 
order 1: U^^^ ~ 0(1) = 1-3). Therefore we can naturally get large 
mixing angles for both atmospheric and solar neutrinos. In our numerical 
analysis, we impose the same constraints for atmospheric, solar neutrino 
solutions and the CHOOZ experiment as in Model I. Then we search for 
the solutions. In Fig. ^, we show their distribution. As we expected, the 
SMA solution is hardly obtained, and it is difficult to make a large enough 
mass hierarchy for the LOW solution. Thus Model II strongly prefers the 
LMA solution. 0.01% of the data sets passed the conditions for the LMA 
solution. The probability to realize the SMA and LOW solutions is very 



small. Therefore a near-future experiment, KamLAND |[T^, will be able to 
test this class of models. In Model II, the rate to realize the LMA solution 
is much smaller than in Model I. The main reason is that the limit on the 
CHOOZ experiment severely constrains the models, since the value of U^^^ 
tends to become of order 1. In Fig. ^ we show the predicted values of U^^^ 
for Model II. As can be seen from Fig. H, the predicted values are so large 
that the future measurement of U^^^ can observe them [Q. Therefore, 
again the measurement of U^^^ is very important. 

Another interesting possibility within Model II is that the inverted hier- 
archical neutrino masses, that is m^g -C < with Sml^^ = ml^ — 
and Sm^^i = ~" ^li^ could be reahzed because of the degeneracy of the 
neutrino mass matrix elements. However, we checked that the rate to realize 
such a possibility is very small; only 2 x 10~^ % of the data sets passed the 
criterion of the inverted hierarchy and the constraints for the LMA solutions. 
This possibility can be tested by future neutrino experiments such as the 
neutrino factory [0, if this could be reahzed. 

Here we mainly stressed the importance of the measurement of U^^^ . 
In addition to this, CP violation and 2f]0i> decay in neutrino physics are 
important to see the distinct features for Models I and II, as pointed out in 
Ref. §. 

In the next section, using the models that satisfied the neutrino con- 
straints here, we will discuss in detail the LFV in the charged lepton sector. 
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Figure 3: Distribution of the points of neutrino solutions in Model II. All 
points satisfy the constraints for atmospheric and LMA solutions in Eqs. ([lID 
and 
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4 Lepton-Flavor Violation 



In the presence of non-zero neutrino Yukawa couplings, we can expect LFV 
phenomena in the charged-lepton sector. Within the framework of SUSY 
models, flavor violation in neutrino Yukawa couplings induces LFV in slepton 
masses even if we assume the universal scalar mass for all scalars at GUT 
scale [|^, |]. In the present models, the LFV is generated in left-handed 
slepton masses since right-handed neutrinos couple to the left-handed lepton 
multiplets. A RG equation for the left-handed doublet slepton masses 
can be written as 



^ dfi \ (i/i 



MSSM 



167r2 



mlflf. + flUml + 2UlmlU + rhljlU + 4^.)! . . , (19) 



where m| and fn\j^ are soft SUSY-breaking masses for right-handed sneu- 
trinos (z>) and doublet Higgs (i^w), respectively. Here {fid{mY)ij / dfi)ussM 
denotes the RG equation in case of the minimal SUSY SM (MSSM), and 
the terms explicitly written are additional contributions in the presence of 
the neutrino Yukawa couplings. In a basis where the charged-lepton Yukawa 
couplings are diagonal, the term (//(i(m|)ij /cZ/i) mssm does not provide any 
flavor violations. Therefore the only source of LFV comes from the additional 
terms. In our analysis, we numerically solve the RG equations, and then 
calculate the event rates for the LFV processes by using the complete formula 
in Ref. |^. Here, in order to obtain an approximate estimation for the LFV 
masses, let us consider one-iteration approximate solution to the LFV mass 
terms {i 7^ j): 

(^,,„|)„^.(£1|K ,og^. (20) 

Here we assume a universal scalar mass (mg) for all scalars and a universal A- 
term {Af = Oomo//) at the GUT scale {Mq = 2 x lO^^ GeV). From Eqs. (D 
and flSI), the solution can be written as 



(6 + flo)'^0 r rDirac rr-Dir-ac* I f |2 i„„ ^( 



{Ami)., ^ - ^g^r" ^"^"^^^"^*!/'-^!' ig- (21) 

Note that large neutrino Yukawa couplings and large lepton mixings U^^^°''^ 
generate the large LFV in the left-handed slepton masses. 
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The component (Am|)32 ((Am|)2i) generates the t ^ fi'y {fi 67) 
process through diagrams mediated by sleptons. The decay rates can be 
approximated as follows: 

(22) 

Here F is a function of masses and mixings for SUSY particles. If the non- 
degeneracy of slepton masses is very small, the function F is approximately 
process- independent. 

Note that the mixing matrix that is relevant to the LFV masses is JJ^^^"-'^^ 
not the neutrino mixing matrix Jj'^^^^ _ The mixing matrix f/^*'"'*^ depends 
on the FN charges of right-handed neutrinos (a, b, c) and 6. In the next 
subsections, we discuss the LFV in some interesting cases. 



r(e, ej7) 



1 f^TT 



(Ami) 



4.1 Model I {6=1) with (a, b, c) = (0, 1, 2) 



First we discuss Model I (5 = 1) with (a, 6, c) = (0,1,2). In this case, the 
neutrino Dirac mass matrix in Eq. (HI) is given by 



( Cse^ B^e' A^e" 
C2e^ B2e Aae 



(23) 



When we diagonalize the Dirac mass matrix, we get the following approxi- 
mate expression for f/^*'''^^. 



U 



Dirac 



I 



\ 



1 



eX 



AiB2-B^A2 

___€B2Cx_ 

A1B2-B1A2 



v/i+|5iMi|2 
1 

y/l+\_Bi/Ai\2 

Bi/Ai 



v/i+|5iMip 

1 



(24) 



v/l+|BiMi|2 y/l+\Bi/Ai\2 j 



where X = C*(y4iA* + BiB;)/Ai{AIB; - A*5*). This expression is a 
leading order in terms of e. An important point is that the lepton mixing 
matrix JJ^'^^"''^ in Eq. (p^ ) also contains a large mixing because of the lopsided 
structure of the Dirac neutrino masses. 

The component (Am|)32 induces the r ^ /i7 process. Taking into 
account the hierarchical structure of neutrino Yukawa couplings {f^i : f^2 '■ 
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fu3 ~ : e : 1), a dominant term of the matrix element (Am|)32 can be 
written as 



(Am 



2^ 



(6 + al)m, 
16^ 



TrDiracTj-Dirac*} f |2 
'^33 '^23 



l/.3plog 



Mr- 



Ad, 



R 



(25) 



The component (Am|)2i, on the other hand, generates the fi —>■ ej process. 
Since U^^'^"''^ is of order e, a leading term of (Am|)2i can be written as 



(Am^ 



(6 



2\ 2 

'^0 / "^q jjDirac jjDirac* 



167r2 



'13 



Mr 



M 



R 



(26) 



The branching ratio for the r — /i7 (/x — 67) process is proportional to 
the square of (Am|)32 ((Am|)2i), as given in Eq. (p^). Thus the important 
parameters for the branching ratios are the matrix elements U^^'^"''^, Uf^^^"''^ 
and U^''^"''^, and the third-generation Yukawa coupling fu^. In Model I, the 
element U^""''^ is expected to be of order 1 because of the lopsided structure 
in the Dirac mass terms as expressed in Eq. (0). In Fig. ^ we show the 



numerical result for U^^^"''^ 

TT Dirac 
'^13 



compared to U^^^ . As expected from Eq. 



is of the order of e, as shown in Fig. |^. The order-1 U^"'""^ and 
the order-e (non-zero) U(^"'"'^ can induce a significantly large event rate for 
/i — >• 67 process if the neutrino Yukawa coupling f^^ is large, as pointed out 
in Ref. 0. 

A magnitude of f^^ depends on the FN charge r. Even at this stage, 
however, we can estimate the ratio of the branching ratios of the /i — 67 and 
r ^ /i7 processes as Br(yU e7)/Br(r — /i7), which is r-independent. If we 
consider a value of T{ei ej7)/r(ej — > CjUiUj), the dependence of the initial 
lepton mass in the process is cancelled. Thus we can expect the following 
relation from Eqs. (^2]) , (p5|), and (|^): 



r(/i- 


e7)/r(/i ev^Ve) 


(Am|)2i 


2 


JjDirac 
•-^13 


2 




r(r- 


^ /i7)/r(^ l^^T^t?) 


(Am|)32 




T T Dirac 
^33 





.(27) 



Taking into account Br(/i ^ ez/^z/g) ~ 100% and Br(r ^iVr^p) — 17%, we 
obtain 



Br(/i — i> 67) 

Br(r ^ /i7) ^ 0l7 



^1 



0.03 




(28) 



This relation is a r-independent prediction of Model I with (a, 6, c) = (0, 1, 2). 
Since the current experimental limits on these processes are Br(yU 67) < 
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Figure 5: U^^^ versus t/^^™^ in Model I with (a, 6, c) = (0, 1, 2). 
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1.2 X 10~^^ and Br(r fi'-f) < 1.1 x 10~^ the process /i — 67 is much 
more sensitive to this class of models. 

In order to discuss the absolute values of the branching ratios for LFV 
processes, we need to fix the FN charge parameter r. In SUSY models, there 
are two interesting cases. In the models with r = 0, it is suggested that all 
the third-generation Yukawa couplings are of order 1. In SUSY models, this 
is realized in the case with the large tan/3 (tan/3 ^ 50). In the models with 
r = 1, all the third-generation Yukawa couplings are of order e, except the 
top Yukawa coupling, which is expected to be of order 1. This is likely to 
be the case with small tan (3 (tan /3 ~ 5) in SUSY models. Therefore, in the 
next subsections, we will consider the two cases. Case 1, r = and Case 2, 
r = 1, to see how large the predicted branching ratios for LFV processes can 
be. 



4.1.1 Case 1, r = 

In Case 1, the models suggest that all third-generation Yukawa couplings at 
the GUT scale are of order 1. Therefore, a large tan/5 is preferred. Here 
we take tan/3 = 50. In our numerical analysis, we impose = mtop at 
GUT scale in the Dirac neutrino mass Eq. (P^D with r = 0, and then we 
require the neutrino mass squared difference to be ^mgg = 3 x 10^^ eV^ 
in order to fix the right-handed neutrino mass scale Mr. We have checked 
that all third-generation Yukawa couplings are of the same order, which is 
about 0.5-0.8. Since the third-generation neutrino Yukawa coupling is as 
large as the top Yukawa coupling, so the large LFV masses are generated 
through the RGE running. Furthermore, since tan/? is very large, the low- 
energy amplitudes of LFV processes are enhanced. Therefore, we can expect 
very large branching ratios for LFV processes fl^. In Fig. ^ we show the 
numerical result for Br(r fi'-f) and Br(/i —>■ ej). Here the circle points 
satisfy the LMA solution for solar neutrinos, the diamond-shaped points the 
SMA solution, and the square points the LOW solution. As can be seen 
from the figure, the branching ratios for the LFV processes do not depend 
on the solar neutrino solutions since the important parameters for n ^ e'j 
and T ^ fi'j branching ratios are only U^^^"''^, U^^'^"''^ and U^^^"''^, which are 
not affected very much by the constraints on solar neutrino solutions, as 
shown in Fig. ^. We also see that the estimation Eq. (^) is approximately 
correct, so that the limit from the /i — 67 process can provide much stronger 
constraints on the models. In Fig. |^, we select one point from Fig. ^ and show 
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Figure 6: Br(/i — > 67) versus Br(r — > /i7) in Model I with (a, 6, c, r) = 
(0,1,2,0). Here we take the left-handed slepton mass to be 300 GeV, the 
Wino mass to be 150 GeV, and e = 0.07. 

-7 tan B= 50 M2 = 150 GeV ag = li>0 

10 \ ^ ^ ^ \ 



p 




200 400 600 800 1000 
left-handed selectron mass (GeV) 



Figure 7: Br(yU —>■ ej) as a function of the left-handed selectron mass in 
Model I with (a, b, c) = (0, 1, 2, 0). 
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the slepton mass dependence of the branching ratio for the fi ^ process. 
In this case, these branching ratios are too large to be consistent with the 
present experimental bounds. Therefore the large region of parameter space 
is already excluded. 

These results illustrate the significant potential of the LFV searches to 
probe the realistic neutrino models. 



4.1.2 Case 2, r = 1 

In models with r = 1, the top Yukawa coupling is expected to be of order 
1, while the other third-generation Yukawa couplings (bottom, tau and tau- 
neutrino Yukawa couplings) are of order e. It is thus likely to be the case with 
relatively small tan/5. Here we take tan/3 = 5, and we impose = mtop 
at GUT scale in the Dirac neutrino mass matrix in Eq. (^) with r = 1, in 
other words, (/jy)33 = eftopAi. We checked that in the case with tan/5 = 5, 
the bottom and tau Yukawa couplings approximately satisfy the condition 
/fe ~ ~ eftop- Q Therefore, as compared with the case with r = 0, the LFV 
masses in Eqs. (pSD and (|^) are suppressed by e^, and then the branching 
ratios for the r — > /i7 and /i — 67 processes are reduced by e^. The numerical 
results for the branching ratios are shown in Fig. |^. Note that the predicted 
branching ratios in all solar neutrino solutions are almost the same, and that 
the relation in Eq. ( ]28|) is approximately satisfied. Thus the — > 67 search 
is much more sensitive to this class of models, unless the future sensitivity of 
Br(r fi'-f) can reach values much below 10~^°. In Fig. |^ we also show the 
branching ratio for yU — > 67 as a function of the left-handed selectron mass. 
As can be seen from Figs. ^ and |^, it is very interesting that the predicted 
branching ratios for /i ^ 67 can be just below the present experimental 
bound (Br(/i — >■ 67) < 1.2 x 10~^^ |T^) in a wide range of parameter space. 



We also get the ratios between branching ratios Br(/i eee)/Br(/i 67) 
and R(/i — e in Ti (Al))/Br(/x 67) as follows: 



Br(yU —>■ 67) 
R(/i ^ e in Ti (Al)) 



6 X 10"^ (29) 
5 (3) X 10-\ (30) 



Br(yU 67) 

In the present models, these ratios are quite predictive since the on-shell 
photon penguin diagram, which induces the /i ^ 67 process, dominates over 

^The tan /3-dependence of the branching ratios is approximately tan^ 
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Figure 8: Br(yU — > 67) versus Br(r yU7) in Model I with (a, 6, c, r) = 
(0,1,2,1). Here we take the left-handed slepton mass to be 300 GeV, the 
Wino mass to be 150 GeV, and e = 0.07. 
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Figure 9: Br(/i — 67) as a function of the left-handed selectron mass in 
Model I with (a, b, c, r) = (0, 1, 2, 1). 
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the other contributions to the /i eee and jj, ^ e conversion processes. 
Therefore the future improvements of the branching ratios, Br(/i 67) ~ 
10-1^ [0, and R(/i ^ e in Al) ~ lO'^^ |6[ (and possibly R(/i ^ e in 
Ti) ~ 10~^^ |T^) will provide a significant impact to this class of models. 



4.2 Model I with (a, 6, c) = (0, 0, 0) 

In this case, the Dirac neutrino mass matrix in Eq. 



is given by 



B2 A2 



(31) 



Then we obtain the mixing matrix JJ^''™'^ as follows: 



U 



Dirac 



( 1 e 

en - 



X- 



X4 



Z 



N+ 
Z 



(32) 



N+ I 



in the leading order of e. Here the coefficients X±, l^±(o), ^(o); and A^± are 
expressed by 



(A,i?;)(A*Q)-(A,4*)(i?|Q) 



(4*i?,)(i?;c,)-(4*c,)(i?;i?,) 
(4*A,)(5|5,)-(A,5*)(A*s,)' 



(4*A,)(5;5,) - (A,i?r )(A*i?,)' 

-(A* A,) + (5*5,) ± J{(A*A,) - (5*5,)}' + 4(A*5.)(A,5;) 



(4*5,), X± = -r±Fo*-^±^o, Ar± = |F±|2 



(33) 



Here summation over i and j (z, j = 1, 2) is assumed: for example, (A*B^ = 
J2i=i-2 ^i^i- The main difference from that in the previous case (see Eq. (^4])) 
is that the matrix elements U^^^"''^ depend on more parameters, since the 



matrix elements (rriiyD 



1, 2, j = 1 — 3) are larger than those in Eq. ( p3| 



However, the important point is the same, that is the matrix elements t/,^"'"'^ 



i,j = 2, 3) in Eq. (|3^) are of order 1 because of the lopsided structure of the 



Dirac neutrino mass matrix, and the element U{^^^°''^ is of order e. 

In the present case, the feature of the LEV is very similar to the previous 
case: Model I with (a, 6, c) = (0,1,2). The hierarchy of neutrino Yukawa 
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couplings is f^i : f^2 '■ fu3 ~ e : e — 1 : 1, which depends on the solar neutrino 
solutions, and hence also the contributions induced by f^2 will be important. 
The leading contribution to the LFV mass (Am|)32 is given by 

(Am|)32 ^ 



In the SMA and LOW solutions, jt/gf ^^t/af ~ 0.5 because of the lopsided 
structure of the Dirac neutrino mass matrix, and (|/i/3p — |/!/2n/(|/i/iP + 
|/!/2p+|/!/3p) ~ 1 since the mass scale for solar neutrinos is much smaller than 
the mass scale for atmospheric neutrinos (i.e. Ifusl \fu2\, \fui\)- On the 
other hand, in the LMA solution, \Ug''^''^Ug'''^''\ and (l/^sp- |^2p)/(|/i.iP + 
|/i/2p + |/i/3p) can be widely spread since other matrix elements of JJ^^^"''^ 
may have a large mixing and the mass scale for solar neutrinos can be close 
to the scale for atmospheric neutrinos. In Fig. [1^, for example, we show the 
distribution of the values for f/^*''"'^ compared with the values for f/^^^'^. In 
the SMA and LOW solutions, U^'"^"-" ~ on the other hand, the values 

for U^"'""^ in the LMA solution are widely distributed. Therefore, the values 
of the branching ratio for r — > /i7 in the case of the LMA solution can be 
much more broadly distributed than in the other cases. 

The dominant contribution to the LFV mass (Am|)2i can be written as 

(Am|)2i ~ 



In the cases of the SMA and LOW solutions, the first term in Eq. (|35D 
is dominant since Ifusl ^ \fu2\ and U^^^"''^ is non-negligible, as shown in 
Fig. 0. In the case of the LMA solution, the second term can be large since 
fu3 ~ fu2 in order to get mass scale for the LMA solution, and U^x^"-'^ can be 
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(6 + a^ 



167r2 



T j-DiracT j-Dirac* \ r i 
^^33 ^23 \ J 1^31 



I T rDirac j rDirac* 



'22 



I/. 



p2\ 



(6 + ag 
167r2 



T rDiracT rDirac* 



23 



1/3 1 ~ \fv2 



log^^- (34) 



M, 



R 



(6 + al)ml r 



167r2 



T rDiraCT rDirac* \ f 1 2 
f^23 '^13 l/i'Sl 



'13 



rDirac j rDirac* 



,T ruiracj T 
+ '^22 ^12 

(6 + al)ml 



\U2? log 



M, 



G 



Mr 



167r2 



T rDirac j rDirac* 
^23 ^13 



(|/.3P 

TrDiracTrDirac*\ f \2 i„_ 



l/.2|- 
Mp 



(35) 



large because of the large mixing for the LMA solution. The matrix elements 
^jjDtracjji^irac^ and | f/|J*^"^ f/j^*'""'^ | are broadly distributed (for example, see in 
Fig. |TOD, and hence the value of the branching ratio for fi —>■ e'y can be widely 
distributed in all cases. 

When |/jy3p ^ \fu2\'^, a ratio of branching ratios can be approximately as 
follows: 

Br(/i — * 67) 1 
Br(r ^ /i7) 0.17 

Taking into account the present experimental limits and the future expecta- 
tion on these processes, the search for /x — >■ 67 can be much more sensitive 
to this class of models than the search for r fi'y. 

In the following subsections, in order to obtain the magnitude of the 
branching ratios for the LFV processes, we consider two cases again: Case 1, 
r = and Case 2, r = 1. 



{Ami 



121 



32 



0.01 



TjDirac 
•-^13 



/0.03 



(36) 



4.2.1 Case 1, r = 

As in the previous models in section 4.1.1, we take tan/3 = 50 in order 
for all the third-generation Yukawa couplings to be of order 1. We assume 



that m3 = rritop at the GUT scale in Eq. (0) with r = 0. In Fig. |TI 
we show the result of the branching ratios for the t fi'j and /i — > 67 
processes. Because of the large neutrino Yukawa coupling, the branching 
ratios of the fi ^ e'-f process are so large that most of the parameter region 
can be excluded by the present limit on /i — >■ 67. As in the models with the 
FN charges (a, 6, c) = (0,1,2), the models with r = 0, which have a mild 
Yukawa unification between the top and tau-neutrino, are almost excluded 
by the bound on /i — > 67 [p!0[] . 



4.2.2 Case 2, r = 1 

In the models with r = 1, as discussed in section 4.1.2, the following relations 
among the third generation Yukawa couplings are implied: ~ /t ~ fu3 ~ 
eftop- Thus we assume that in Eq. ( PT| ) with r = 1, 1713 = rritop at the GUT 
scale and take tan/? to be 5 in order to realize the relation. As compared 
to the case with r = 0, the branching ratios for r fi'-f and /i — > 67 are 
suppressed by e^, with the additional suppression due to small tan/3. In 
Fig. |12| we show the result of the branching ratios for r — > yU7 and fi — > ej. 



22 



Dirac 
U23 



1 

0.8 
0.6 
0.4 
0.2 








J -.Dirac 
Ui3 -1 

10 



10 



10 



•LMA 
♦ SMA 
■LOW 



10 



0.2 



• ♦ 
% 



S 

♦ • •• ^tl^ 



0.4 0.6 

MNS 



0.8 



10" 



10 

MNS 

Ue3 



10 



• LMA 






♦ SMA 






■ LOW 






■ 




m 




♦ • 

■ ♦ 


•♦ 



Figure 10: U^'^^^ versus f/^^™^ in Model I with (a, b, c) = (0, 0, 
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Figure 11: Br(/i —>■ 67) versus Br(r fi'-f) in Model I with {a,b,c,T) = 
(0,0,0,0). Here we take the left-handed slepton mass to be 300 GeV, the 
Wino mass to be 150 GeV, and e = 0.07. 
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Figure 12: Same as Fig.|Tl] except for Model I with (a, b, c, r) = (0, 0, 0, 1). 
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For the jj, ^ 3e and ji ^ e conversion processes, we have checked that 
the relations in Eqs. (pQ]) and (|30D are held. Therefore, interestingly, the 
predicted event rates for /i — > 67 and ^ e conversion can be as large as 
those the future experiments can reach. On the other hand, for r — >■ /i7 
process, the significant improvement of the branching ratio sensitivity of 
much below 10"^'^ will be needed in order to reach the predicted values. 



4.3 Model II {5 = 0) with (a, b, c, r) = (0,1, 2, 1) 

Let us discuss the LFV in Model II (5 = 0). First we consider the case with 
hierarchical right-handed neutrino masses, that is the case with (a, b, c) = 
(0, 1, 2). The Dirac neutrino mass matrix is given by 

m^D = m^e C^e B^e Age . (37) 
V C, B, A, J 

Since the matrix elements {muD)3i/{^3^) {i = 1^3) are of order 1, we can 
expect that all elements of the mixing matrix JJ^"'"-^ would be of order one: 
jjDirac ^ 0(1), even if we imposed the constraints of neutrino parameters 
on the neutrino mixing matrix JJ^^^^ _ This structure of the mixing matrix 
jjDtrac jg q^j^g distinct from that in Model I, as seen in Eqs. ( P^ and (P^. 
For example, in Fig. 0, we show U{^^^°''^ compared with f/^*''"^. Even though 
the CHOOZ limit \U^^^^\ < 0.15 was imposed, the component t/j^*™'^ can 
be of order 1. This large 11^^^°''^ is a very important feature for the LFV in 
the present models. 

Because of the FN charges of right-handed neutrinos, the structure of 
neutrino Yukawa couplings is hierarchical {f^i : fu2 : /1/3 ~ : e : 1)- 
Therefore, only the third-generation Yukawa coupling f^^ is important for 
contributions to LFV masses. The dominant terms of the matrix elements 
((Am?)32) and ((Am^)2i) are given by 



(Am|)32 ^ -^^^T^Ui--Ui"-nU'log^, (38) 



R 



{AmlU - - \^fj"' Ui--UR-nU' log (39) 
Then a ratio of branching ratios Br(/i e7)/Br(r /I'f) are approximately 
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Figure 13: t/^f versus Ug'"'"^ in Model II with (a, 6, c) = (0, 1, 2). 



expressed as follows: 

Br(/i — >• 67) 



Br(r ^ /i7) 0.17 



(Am' 
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(Am?) 



32 



TT Dirac 
^13 



T T Dirac 
'^33 



(40) 



Since |f/{^*''°'^| ~ |f^33*'^"'^|; as shown in Fig. O, the branching ratio for /i — > 67 
can be as large as or even larger than that for r — *• /X7. As discussed in Model 
I, the models with r = 0, in which a mild Yukawa unification between the 
top and tau-neutrino is realized, are almost excluded by the present /i ^ 67 
bound 



10 



Therefore, in this section, we only present the result in the 
models with r = 1. 

In Fig. |l^, we show the numerical result for these branching ratios. The 
predicted branching ratios for — >■ 67 are larger than those in Model I, since 
the mixing element Ug^^"''^ is much larger than that in the Model I. On the 
other hand, the branching ratios for t —>■ fij are almost the same as those 
in Model I. f\ As expected, the predicted branching ratios for fi e'-y are as 
large as those for r fij. Note that for the fi —>■ 3e and fi ^ e conversion 



^If the inverted neutrino mass hierarchy is reahzed, as discussed in section 3.2, neutrino 
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Figure 14: Br(/i — 67) versus Br(r ^ /i7) in Model II with (a, 6, c, r) = 
(0,1,2,1). Here we take the left-handed slepton mass to be 300 GeV, the 
Wino mass to be 150 GeV, and e = 0.07. 



processes, the relations in Eqs. (|^) and ( ^OD are predicted. Therefore the 
search for LFV in muon processes is much more sensitive to this class of 
models too. Even at present, some of the points in Fig. |l^ are already 
excluded by the current experimental limit on /x — 67 [Br(yU — > 67) < 
1.2 X 10~^^ [mi]- The future proposed improvement of the limit on LFV 
in muon processes will definitely provide a significant test of this class of 
models. 



4.4 Model II {6 = 0) with (a, b, c, r) = (0, 0, 0, 1) 

Next we consider the models in which the heavy right-handed neutrinos are 
nearly degenerate, that is the models with (a, b, c) = (0, 0, 0). Again here we 

Yukawa couplings /i^i and f^2 can be so large that they can contribute to the LFV masses. 
In this case, we have checked that the branching ratios for the LFV processes, especially 
for /i 67, could be more enhanced because of the additional contributions from f^i and 
/i/2 to the LFV masses. 
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Figure 15: Ug''''"' versus Ug"'"" in Model II with (a, 6, c) = (0,0,0). 



only discuss the case with r = 1. 

The Dirac neutrino mass matrix is given by 



C2 B2 A2 
V C, 5i A, 



(41) 



As in the previous case of Model II, we expect Ujj'^^"''^ 
case, if coefficients Aj, Bj and (7,- are all real, JJ^^"^""^ = 



~ 0(1). In this 
UMNS_ '^Yms, as 

shown in Fig. 13, the component t/^*''"'^ tends to be smaller than that in the 
previous case, because the constraint \U^'^^\ < 0.15 slightly affects the value 
of Ug^^"''^. However, all components in the Dirac mass matrix Eq. ( p] ) are of 
order one, and all Yukawa couplings f^i {i = 1-3) can contribute to the LFV 
in slepton masses. Therefore, branching ratios of r ^ /i7 and fi e'-f can be 
larger than those in the previous case of Model II. We present the numerical 
result in Fig. 0. As can be seen in Fig. 



m 



the branching ratio for /i — > 67 
because all components V^^'^"''^ can be 



can be as large as that of r 
of the same order. Many of the points are already excluded by the present 
experimental limit on the /i — > 67 process, and the future improvement will 
be able to test this class of models. 
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Figure 16: Br(// — > 67) versus Br(T jji'y) in Model II with (a, 6, c, r) = 
(0,0,0,1). Here we take the left-handed slepton mass to be 300 GeV, the 
Wino mass to be 150 GeV, and e — 0.07. 

5 Conclusions 

We discussed the neutrino masses and LFV in SUSY models with lopsided 
FN charges. For neutrino physics, we stressed the importance of the mea- 
surement of U^^^^^ in both Models 1 and 11. In Model 11, the LMA solution is 
strongly preferred so that the near-future KamLAND experiment could play 
a significant role for this class of models. 

We also analyzed the LFV in detail within this framework. The present 
experimental limits on LFV processes almost excluded the models with r = 
in which a mild Yukawa coupling unification between the top and tau- 
neutrino is realized at the GUT scale. In the models with r = 1, in which 
the tau-neutrino Yukawa coupling is suppressed by e compared with the top 
Yukawa coupling, the predicted branching ratios for r — > //7 are as large as or 
much less than 10~^°, and hence a significant experimental improvement of 
the limit on the branching ratio is needed in order to reach the predictions of 
this framework. On the other hand, the predicted branching ratios for 
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67 and /i — > e conversion in nuclei can be as large as those the future proposed 
experiments can reach. Therefore, the future experiments can provide the 
significant test of this realistic framework for neutrino masses as well as 
charged lepton and quark masses. 
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